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Self-similar spectrum in effective time independent Hamiltonians for kicked systems 

Rashmi Jangid Sharma,^’0 Jayendra N. Bandyopadhyay,^’[i| and Tapomoy Guha Sarkar^’d 
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We study multifractal properties in the spectrum of effective time-independent Hamiltonians 
obtained using a perturbative method for a class of delta-kicked systems. The evolution operator 
in the time-dependent problem is factorized into an initial kick, an evolution dictated by a time- 
independent Hamiltonian, and a final kick. We have used the double kicked SU (2) system and the 
kicked Harper model to study butterfly spectrum in the corresponding effective Hamiltonians. We 
have obtained a generic class of SU{2) Hamiltonians showing self-similar spectrum. The statistics 
of the generalized fractal dimension is studied for a quantitative characterization of the spectra. 
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Hamiltonian systems undergoing periodic delta kicks 
are studied extensively as a generic model for classi¬ 
cal and quantum chaos [I|. This has found new rele¬ 
vance in the possibilities of engineering such systems us¬ 
ing ultra-cold atoms 0. In the traditional approach, 
such time-dependent systems are theoretically studied 
using the Floquet analysis whereby the quasienergy spec¬ 
trum is investigated for the signature of quantum chaos 
[l[ and quantum criticality with varying parameters of 
the Hamiltonian ii- Quantum chaos studies have also 
shown the existence of fractal butterfly patterns in the 
quasienergy spectrum of periodically driven systems [1-Q 
indicating an infinite number of quantum phase transi¬ 
tions 0. These systems are particularly interesting in 
the fact that though their classical phase space dynamics 
may be chaotic, the quantum quasienergy spectrum does 
not follow the celebrated Bohigas-Giannoni-Schmit con¬ 
jecture 0- We focus on two such systems. The first one 
consisting of driven SU (2) operators is also known as the 
double kicked top model. This system is of interest as it 
finds realization through driven two-mode BEG systems 
d, 0 ■ Secondly, we look at the kicked Harper model. 
The time-independent version of the Harper model rep¬ 
resents the behavior of electrons in periodic lattice in the 
presence of a constant magnetic field in the tight-binding 
nearest neighbor approximation [ll[ll. 

It is possible to construct a time-independent effective 
approximate Hamiltonian for such time-dependent sys¬ 
tems when the frequency of the periodic driving is large. 
Traditionally, the effective Hamiltonian is obtained from 
the Floquet operator using the Cambell-Baker-Hausdorff 
(GBH) or Trotter expansion. It has been shown that the 
GBH method to study the kicked systems suffer intrin¬ 
sic flaws and an alternative formulation [l^, [l^ is better 
suited for more accurate analysis of such systems . The 
effective Hamiltonian thus obtained is found to mimic the 
exact time-evolution for a large range of parameter val¬ 
ues. 
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In this paper, we have investigated for fractal spec¬ 
trum in the effective time-independent Hamiltonian ob¬ 
tained from the Floquet operator using the perturbative 
method used in earlier works H] • We have used the 
double kicked SU{2) system [y,[3 and the kicked Harper 
model Q to study fractal spectrum in the corresponding 
effective Hamiltonians. We have obtained a generic class 
of SU (2) Hamiltonians showing self-similar spectrum in 
finite dimension. These models revealing butterfly spec¬ 
tra also contain the Harper model as a special case. We 
study in detail the (multi)fractal properties of the eigen¬ 
values and eigenstates of the effective Hamiltonians for all 
the models considered. The statistics of the generalized 
fractal dimension is studied to quantitatively understand 
their scaling behavior. 

A general time-dependent problem where H(t) = Hq -|- 
V{t), with a time-periodic potential V{t) = V{t + T) of 
periodicity T has a Floquet operator JFit) which cor¬ 
responds to the time-evolution operator for one time- 
period. The traditional approach to extract an effective 
static Hamiltonian one writes T = exp(—liJeffT) and 
uses the GBH expansion to read out iJeff up to any order 
in T ~ 1/lo. This method, however has been found to 
suffer from several inadequacies. 

The method used in Refs. [3 [II expresses the time- 
evolution operator U(ti —^ t/) between time instants ti 
and tf = ti + T, as a sequence of operations consisting 
of an initial kick followed by an evolution under a time- 
independent Hamiltonian and a final ‘micro-motion’ 

UiU^tf)=Witf)e-^^‘^<^UiU) ( 1 ) 

where U(t) = such that F{t) = F{t + T) with 

zero average over one time period. For high-frequency 
pulsing, the operators i/eff and Fit) can be expanded as 
a perturbation series in 1/uj of the form 

OO - OO -| 

"•» = E U) = E (2) 

n=0 n—1 

This ansatz along with Eq. m can be used to obtain 
iJeff and F{t) up to any desired accuracy. In this method. 
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can also be obtained from a different driven SU (2) Hamil¬ 
tonian of the form 


+ 00 

'H = 'Ho + V ^ S{t — nT), where V = aJx 

n— — co 


and 


^0 = exp 




-f- h.c. 


( 5 ) 


This is a single kicked system whose Floquet operator, 
given in Eq. o, matches exactly with that of the double 
kicked top and thereby exhibits interesting Cantor set 
properties in the quasienergies spectrum. The possibility 
of experimental realization of this system has been stud¬ 
ied [l^. We are interested in the spectral properties of 
the effective approximate static Hamiltonian correspond¬ 
ing to this system. The effective Hamiltonian is given by 


FIG. 1. (Color online) (a) Folded energy spectrum of the effec¬ 
tive time-independent Hamiltonian of the double kicked top 
showing butterfly pattern, (b) Multifractal scaling exponent 
showing a linear dependence on q for rj/j = Gr = {V5— l)/ 2 . 
Insets: Self-similar distribution of energy spectrum on differ¬ 
ent scales, (c) Eigenvectors in eigenbasis. D 2 values are 
0.002, 0.218, 0.511, and 0.855 from top to bottom. Localized 
states have smaller D 2 . (d) Corresponding multifractal be¬ 
havior for the same eigenstates. Slopes in the linear region 
are -0.003 (A), -0.314 (x), -0.489 (□), and -0.771 (o). 


the average time-independent part is retained in Hgg and 
all the time-dependence is pushed to the operator F{t) 
at each order of perturbation. The convergence of the 
perturbation series has been surmised in earlier works 
MM- The periodic potential V {t) may be expanded in 

a Fourier series asV{t) = Vo+J2 { J. 

The truncated series for and F{t) up to 0{l/w'^) can 
be expressed in terms of the Fourier coefficients of V{t) 
[14 fl^ . We use this as the general expression for the 
effective Hamiltonian for periodically driven systems. 

We consider the double kicked top model with the 
Hamiltonian in 


- _ 2a - 7? -p 


+00 r 

E 


S { t — nT — 


— S{t — nT) 


The JiS here represent SU{2) generators in the d = {2j + 
1) dimensional Hilbert space. The corresponding Floquet 
operator is given by [1, Q 


F = exp -b h.c. I exp , 

(4) 

where J+ denotes the operator {J^ + iJy)/2. The 
quasiperiodic nature of the factor for ir¬ 

rational values of the parameter rj/ j leads to interesting 
spectral properties The above Floquet operator 


^ ^ V 1 ^ ^ ^ I \ 

H-- + V] ( ^ ^ 

= Ho + ^ +A[[V,Ho],V] 


\n—l 


( 6 ) 


r 24‘ 


Fig. [TKa) shows the folded energy spectrum of the 
Hamiltonian 'HeS as a function of ^ = rj/iTj for a = 1/j 
where we have chosen an even value of spin j = 20. We 
note that odd values of j would bring about changes in 
the spectrum. The spectacular butterfly appearance for 
the static approximate eigenspectrum is in remarkable 
agreement with the quasienergy spectrum of the original 
double kicked top [y, Q- The spectrum shows qualita¬ 
tive similarity with the Hoftstadter butterfly [T^ owing 
to the presence of the quasiperiodic term. This feature 
is however along the off-diagonal nearest neighbor band 
and is therefore different from the Harper/Hofstadter 
case where it appears along the diagonal. In order to 
study multifractality of the energy spectrum, we set 
at an irrational value of the golden ratio Gr = 1)/2 



FIG. 2. (a)-(b): The distribution of D 2 and D 5 for all the 2j + 
1 = 2001 eigenstates, (c) The distribution of the slope of the 
Tq VS. q curves in the linear region, (d) The distribution of the 
participation ratio indicating that the bulk of the eigenstates 
are localized. 
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and a large value of j = 2500. To study the statistical 
property of the energy spectrum, we consider histogram 
of eigenvalues for different scales. The level distribu¬ 
tion p{E) exhibits remarkable self-similarity as seen in 
insets of Fig. [TJb). The quantitative measure of the self¬ 
similar behavior is done by using the generalized fractal 
dimension. Dividing the full range of the energies into 
N bins of size s each we use the standard box-counting 
to obtain the probability Ps{i) of finding a given energy 
eigenvalue in the i-th bin. The scaling exponent Tg is 
related to the q-th moment via the partition function 

N 

^q{^) = S Ps{iY ^ s'^"- The generalized fractal dimen- 
2 — 1 

sion is defined as Dg = Tg/{1 — q)) [I3. Figure [TJb) 
shows that for large values of j the spectrum indeed 
shows multifractal behavior with Tg linearly decreasing 
with increasing q with a slope —0.871. The quantitative 
measure D2 for this spectrum has the value 0.913. To 
look for fractal behavior in the eigenvectors of "H, we con¬ 
sider four eigenvectors shown in Fig. [TJc) with different 
localization properties. The highly localized eigenvectors 
manifest as the ones which have very sharp support over 
a narrow band of the index m labeling the components in 
Jz eigenbasis. To define the scaling exponent fg for the 
n-th eigenvector we consider its components {cm^} and 
divide the total dimension d = 2j + 1 into M partitions, 
and define p{i) = ^ |cm^ P where the summation extends 
over the components m in the i-th partition. The scaling 
exponent Tg is given by ~ where I = d/M. 

i 

Figure [Hd) shows the scaling of Tg with q. The local¬ 
ized states, as is expected, have very feeble multifractal 
behavior. The fractal dimension D2 has higher values 
for states which are more delocalized. These delocalized 
states also exhibit a faster approximately linear fall of 
the scaling exponent Tg with q. 

To study the statistical properties of the eigenstates, 
we consider the distribution of various quantifiers of frac¬ 
tal behavior. Fig. [2l^a)-(b) show the distribution of 
D 2 = —T 2 and D 5 = —Ts/d over all the eigenstates. 
We find that about ~ 40% of the eigenstates have very 
small values of the fractal dimensions. These states do 
not exhibit any fractal nature. Significant fractal behav¬ 
ior (0.1 < {D 2 OTD 5 ) < 0.8) is exhibited by relatively 
small fraction of the eigenstates. Fig. [2Kc) shows the dis¬ 
tribution of the slopes p of the Tgvs.q curves. We find 
that the large number of eigenstates which correspond 
to small values of D2 and also do not show multi- 
fractality and for these states ^ ^ 0. We conclude that 
the eigenstates which exhibit fractal properties are small 
in number and these states also exhibit multifractality. 
The participation ratio (PR) for the n-th eigenstate de¬ 
fined as PR = measures the number of 

basis states over which the given state has significant 
support. Fig. HKd) shows the distribution of the PR 
for all the eigenstates. We see that about 50% of the 
eigenstates have PR less than 20 which means that, for 


these states, out of d = 2001 basis states, about 1980 ba¬ 
sis states do not have any component. These states are 
hence extremely localized. An important feature of the 
eigenstates of the effective Hamiltonian is the existence 
of a dominant proportion of localized states. These are 
also the states which do not exhibit fractal property. 

Having studied the self-similar spectrum of the effec¬ 
tive Hamiltonian in Eq. we propose a general form 
of such Hamiltonians constructed using SU (2) operators. 
We consider a Hamiltonian of the form 


H = a X. 


C cos(A) -I- h.c. 


(7) 


where A = ^ (|TO)(m-|-1|-|-|m-|-1 )(to|) in the Jz 

m=-j 

eigenbasis {|to)} and X = iq{2Jz + 1)/2X We note 
that the operator A is tri-diagonal in this representation 
with diagonal elements zero and off-diagonal elements are 
unity. The presence of the cosine term may lead to frac¬ 
tal spectrum for irrational values of rj/j. Fig. [3] shows 
the energy spectrum for several choices of parameters a, b 
and operator C as summarized in Table HI 
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a 
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-0.800 

3(c) 

a 

0 

P 

-0.756 

3(d) 

0 

a 

OlJx 

-0.833 

3(e) 

a 

ea 

al 

-0.851 

3(f) 

0 

a 

al 

-0.579 


TABLE I. 

All these cases reveal different forms of butterfly spec¬ 
trum. Fig. [SKa) is in fact very similar to the effec¬ 
tive Hamiltonian for the double kicked top given in Eq. 
(©• We note that Eig. EKa), (b) and (d) correspond to 
the Hamiltonian having only non-zero super- and sub¬ 
diagonal elements whereas m c), (e), and (f) have both 
diagonal and off-diagonal elements. Eig. E])!) corre¬ 
sponds to the Harper/Aubry-Andre Hamiltonian where 
the super- and sub-diagonal elements are unit y a nd the 
cosine modulation is along the diagonal [111,[ij,[l3 ■ Mul¬ 
tifractal property of all these energy spectra is shown in 
Fig. Elg)-(h). Table H] summarizes the values of the slope 
p for all the Tg — q curves in Fig. EKh). Though the spec¬ 
tra look visually different from each other their fractal 
properties show remarkable similarity as is quantitatively 
encapsulated in Fig. Elh). 

The Hamiltonian in Eq. o comprising of SU (2) gen¬ 
erators with a = 0, 6 = a and C = al has a close resem¬ 
blance with the Hamiltonian representing non-interacting 
electrons moving in a 2D periodic square lattice with an 
external magnetic field. This variant of the usual Lan¬ 
dau level problem where translation invariance of the lat¬ 
tice is broken in the presence of a constant magnetic field 
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FIG. 3. (Color online) (a)-(f): The butterfly spectrum for the Hamiltonians summarized in TableH) (g): Shows the multifractal 
property of these spectra. 


maps to the celebrated Harper /Aubrv-Andre equation in 
the tight-binding approximation [l^, ll2, Il3| ■ The Hamil¬ 
tonian consists of an uniform nearest-neighbor hopping 
contribution and an onsite potential varying periodically 
with lattice site: 

L 

^ 2 cos(27rnCT) |n)(n| -|- (|n)(n -b 1| -b h.c.) (8) 

n— 1 

where the summation extends over all the lattice sites. 
The corresponding energy spectrum for irrational cr gives 
the Hoftstadter butterfly. If the onsite term is switched 
on and off at regular interval of time T one has the 
kicked Harper model. The Floquet analysis for the 
kicked Harper model has been studied We investi¬ 
gate the effective time-independent approximate Hamil¬ 
tonian obtained using Eq. ®. We have assumed TLq = 

L ^ L 

a(|n)(n-bl| -b h.c.) and V = Y 2q; cos(27rn(T) |n)(n| 

n—1 n—1 

which yields -b where 

^ 1 ^ 

^ cos^(27rncr) (|n)(n-b 1| -b h.c.). (9) 

n—1 

The correction term in Eq. accounts for the kicked 
nature of the onsite potential in the effective static ap¬ 
proximation. 

Eigure IDJa) and (b) show energy spectrum of H^^'> 
and , respectively. Both show almost similar but¬ 
terfly pattern. The self-similar nature of the spectra for 
tr = Gr is evident in Fig. |lKc)-(f) which shows the spec¬ 
tral density at different energy scales. The multifractal 
nature of the spectrum of is shown in Fig. SKg). 
The linear fall of Tq with q has a slope of —0.597. The 

difference of Tq for the spectrum of and as 

shown in the inset indicates that the multifractal prop¬ 
erties of the eigenvalues are not significantly different for 
these Hamiltonians, with |Ar|max ~ 0.02. The fractal 


properties of the eigenvectors however show remarkable 
difference. Figure (Hh) and (i) show the distribution for 
D 2 and Us for the eigenvectors of the two Hamiltonians. 
The mean D values for the eigenvectors are significantly 



FIG. 4. (Golor online) (a) Eigenspectrum of Harper Hamil¬ 
tonian (b) Spectrum of the effective Hamiltonian for 

kicked Harper model ■ (c)-(d) and (e)-(f) The spectral 
density at two different energy-scales revealing self-similarity 
for and , respectively, (g) Multifractal behavior of 

the eigenvalues of ■ Inset shows the deviation from 
(h)-(i) The distribution of D 2 and Ds for the eigenvectors. In 
each of (h) and (i), the flat distribution in the left corresponds 
to and the sharp distribution on the right corresponds 

to . (j) The multifractal property of few typical eigenvec¬ 
tors with D 2 values 0.248 (o), 0.497 (□), and 0.684(A). (k) 
The distribution of slope Jl of the Tq — q curves in the linear 
regime for all the eigenstates. The steeper distribution on the 
right of the figure corresponds to and the one to the left 
corresponds to 
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different. The eigenstates of the show less variance 
around the mean fractal dimension as compared to the 
eigenstates of Figure SJj) shows the multifractal 

behaviors of three typical eigenvectors of . The dis¬ 
tribution of the slopes /I of the Tq — q curves are com¬ 
pared for the eigenstates of and in Fig. Hl^k). 
Hence we find that the effective static Hamiltonian for 
the kicked Harper model though gives a very similar en¬ 
ergy spectrum as the original Harper system, the fractal 
properties of the eigenvectors are considerably different 
for the two cases. 

We conclude by noting that a wide class of kicked 
systems with Floquet butterfly spectrum also show self¬ 
similar behavior in the energy spectrum of their corre¬ 
sponding effective static Hamiltonian. The multifractal- 
ity in the energy eigenstates are also found to be an use¬ 
ful quantifier to distinguish the self-similar properties of 
the spectrum for different Hamiltonians even when their 
eigenvalues show insignihcant difference in fractal prop¬ 
erty. The effective Hamiltonians, though approximate, 
can be used to study statistical properties of such self¬ 
similar spectra for a wide range of time-dependent prob¬ 
lems. 


I. SUPPLEMENTARY MATERIAL: 
SELF-SIMILAR SPECTRUM IN EFFECTIVE 
TIME INDEPENDENT HAMILTONIANS FOR 
KICKED SYSTEMS 

We consider a general time-dependent Hamiltonian 
H{t) = Hq + V (t), with a time-periodic potential V {t) = 
V{t + T) of periodicity T has a Floquet operator J^{t) 


which is the time-evolution operator for one time-period. 
The method used in Refs. [H, [l^ factors the time- 
evolution unitary operator U(ti —> tf) between times ti 
and tf = ti + T, as a sequence consisting of an initial 
kick followed by an evolution under a time-independent 
Hamiltonian and final kick [T3| . 

Uiu ^ tf)=U\tf)e-^^‘^^^U{ti) ( 10 ) 

where U(t) = so that F(t) = F{t + T) with van¬ 

ishing average over one time period. For high-frequency 
forcing, the operators i/eff and F{t) are expanded as a 
perturbation series in 1 /w given by 

CO 1 OO -| 

= E P{t) = ^ -F(“>. (11) 

n=0 n—1 

This along with Eq. m can be used to obtain Hgg 
and F(t) up to any desired order of perturbation. At 
each order of perturbation, the average time-independent 
part, in this method, is retained in Ffeff and all the time- 
dependence pushed to the operator F{t). The conver¬ 
gence of the perturbation series is to be checked on case 
by case basis [H, [l^ . Expanding the periodic potential 
V(t) in a Fourier series we have 

CO 

V{t) = Vo + Y^ . (12) 

n—1 

In terms of the Fourier coefficients, the truncated series 
for HeS and F{t) up to 0{l/uP‘) may be written as 


1 CX) 1 1^1 

= ffo + k'o + - E - ^ E + h-C.) 

n—1 n—1 

. CO . 

3(()2 E —2[l4, [V_m, Wi-n]] + h.C.^ 

n,m—l 

CO . . OO - 

F{t) = —y - (We™""* - + — V — ([W, i?o + f>o] e™""* - h.c. 

n—1 


(13) 


-h 


E 


1 


n—1 


2 iuj‘^ ^ n(n + m) 

n,m—l ^ 


1 1 

n/m=l ' ^ 


This general expression for the approximate effective 
static Hamiltonian for periodically driven systems is used 
in the article. 
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